Diffusion and the Mesoscopic Hydrodynamics of Supercooled Liquids 
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The description of molecular motion by macroscopic hydrodynamics has a long and continuing 
history. The Stokes-Einstein relation between the diffusion coefficient of a solute and the solvent 
viscosity predicted using macroscopic continuum hydrodynamics, is well satisfied for liquids under 
ordinary to high temperature conditions, even for solutes as small as the solvent. Diffusion in 
supercooled liquids near their glass transition temperature has been found to deviate by as much 
as 3 orders of magnitude from that predicted by the Stokes-Einstein Relation M. Based on the 
random first order transition theory j^-Q, supercooled liquids possess a mosaic structure. The 
size- and temperature-dependence of the transport anomalies are quantitatively explained with an 
effective medium hydrodynamics model based on the microscopic theory of this mesoscale, mosaic 
structure. 
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According to macroscopic continuum mechanics and 
Einstein's Brownian motion theory, the diffusion con- 
stant of a probe particle is determined by the viscosity of 
the liquid through which it moves by the Stokes-Einstein 
Relation 
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where Df~ E is the translational diffusion coefficient, r\ 
is the viscosity, and R is the radius of the probe. This 
relation must hold if the solute is much larger than any 
relevant length scale in the liquid. Remarkably it is well 
satisfied even for self diffusion, i.e., where the radius of 
the solvent is the same as the solute. It is actually not a 
bad relation even for diffusion in moderately supercooled 
liquids (T > T g ) (§-[?]]. Viscosity can vary by as much as 
12 orders of magnitude and diffusion coefficients vary by 
nearly as much. Nevertheless there is a noticeable dis- 
crepancy with the macroscopic hydrodynamic prediction 
in the vicinity of glass transition temperature T g . In su- 
percooled liquids around T g , the translational diffusion 
is enhanced by as much as 3 orders of magnitude out 
of the 14 orders of magnitude change of viscosity |@MEJ ! 
Cicerone and Ediger have found the deviation depends 
strongly on the size of the probe and the temperature at 
which the measurement is taken JjJ. For deviations as 
large as a factor of 1,000, small adjustments of hydro- 
dynamic boundary conditions arising from the details of 
the probe-liquid interface will not suffice. It is reason- 
able to conclude there is an additional length scale in 
a supercooled liquid to explain the discrepancy. Such 
a length scale arises naturally in the random first order 
transition theory of glasses. In this paper we will sketch 
the predictions of that theory for the deviations from 
the Stokes-Einstein law expected in supercooled liquids. 
Ediger and others recently have used the size-dependent 
deviations to infer the size of heterogeneous domains in 
supercooled liquids M]. In the published models, glass- 



forming liquids at low temperature are modeled as a two- 
component medium. One component is a solid-like back- 
ground in which the diffusion is slow or impossible; the 
other component is a liquid-like excitation in which stress 
can be released quickly and the probe can move faster. 
By adjusting the parameters like the lifetime, appear- 
ance rate, internal viscosity and the size of the liquid-like 
domains, Stillinger and Hodgdon and Cicerone, 

Wagner, and Ediger |f| have shown it is possible to get 
significant deviations from the Stokes-Einstein law. It 
has been appreciated |L3) that there might not be such 
a clear cut distinction between two components; instead 
there should be a distribution of heterogeneities. A sensi- 
ble question to ask is then how to derive this distribution 
from microscopic theory and thereby to explain quanti- 
tatively the decoupling between diffusion and viscosity in 
supercooled liquids. 




FIG. 1. A probe particle (shaded circle) with radius R is 
shown in a mosaic-like supercooled liquids. At the instant 
shown, the particle is moving inside one of the domains with 
size L. The viscosity inside the domain is denoted as r\i and 
r\ for the averaged outer region. 
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Recently, we developed a microscopic theory of acti- 
vation barriers and fragilities in supercooled molecular 
liquids based on the random first order theory for the 
liquid-glass transition Q). Several remarkable quantita- 
tive regularities emerge from theory. The relation be- 
tween the typical activation barrier in a supercooled liq- 
uid and its configurational entropy is quantitatively ob- 
tained Q] . Furthermore the theory yields the distribution 
of activation barriers explaining the j3 parameter of the 
stretched exponential relaxation functions ]I3 . The the- 
ory also yields the size scale of heterogeneous relaxing 
domains |4|]l4|]. We are now ready to investigate the ef- 
fect of such heterogeneities on molecular hydrodynamics 
in supercooled liquids. To set the stage we briefly review 
the random first order theory and the microscopic theory 
of barriers. 

The Random first order transition framework for the 
liquid-glass transition suggests supercooled liquids relax 
in an activated manner below a certain temperature, Ta- 
-The system has to overcome some free energy barrier 
to reach another metastable state. The activation bar- 
rier arises from a competition between the favorable in- 
crease in configurational entropy (the large number of 
other minimum energy configurations the local region 
could freeze into) and an unfavorable interfacial energy 
(the energy cost of forming another state in the original) . 
This interfacial energy reflects the cost of a domain wall- 
like excitation where atoms are strained and are in far 
from local energy minimum configurations. Mathemat- 
ically, the free energy change in forming a droplet with 
radius r is 
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where s c is the configurational entropy density driving 
the random first order transition. The surface tension is 
strongly renormalized by the multiplicity of states which 
can wet the interface, much like in random field Ising 
model, fi(r) = coIt 1 ) 1 ^ 2 @> wnere r o is the interparti- 
cle distance. The short range surface tension cr is nearly 
universal since it depends on the degree of vibrational lo- 
calization in aperiodic minima. This is usually expressed 
as a universal Lindemann-like criterion for vitrification 
The free energy barrier determined by Eq.(2) can be 
expressed as 



37ro-gr 
Tsc 



(3) 



When the Lindemann value is substituted into the mi- 
croscopic expression for cto this results in an excellent 
description of the activation barriers in a wide range of 
liquids. The size of the droplet (or domain) is determined 
by the condition F(L) — giving 
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This suggests regions of size L are separated by domain 
walls and can reconfigure nearly independently. The con- 
figurational entropy, the driving force for reconfigura- 
tions, fluctuates around a mean value = Ac p (T) T ^ K , 
where Acp(T) is the specific heat jump per unit volume 
at the transition, with a magnitude Ss c = y^CijAc^/W, 
where V$ = §7T-L 3 is the volume of a typical domain. This 
results in a corresponding variation in free energy barri- 
ers and size for each mosaic element. This therefore gives 
a distribution of relaxation times and domain sizes. This 
leads to a quantitative relation between the Kohlrausch- 
William- Watts exponent (3 and the liquid fragility jTJJ] . 

The intricate mosaic structure of the supercooled liq- 
uid state presents a challenging problem in hydrodynam- 
ics. We can think of each region as having its own vis- 
cosity r]i determined mostly by the local relaxation time 
Ti. The distribution of Ti is determined by the entropy 
fluctuations. Of course there will be fluctuations in the 
higher frequency elastic modulus of each region, yi, but 
these should be small. Thus r\i «< y > t,;. If the probe 
is as large or larger than a typical domain it will be sur- 
rounded several cells of the mosaic with varying rji. But 
the effect of this will be essentially the same as a homo- 
geneous medium with the average viscosity rj —< rji >. 
Only probes smaller than L will be affected. The situa- 
tion should be well described by a two-zone fluctuating 
viscosity hydrodynamics: a viscosity rji locally, which has 
a fluctuating value embedded in an infinite region with 
the average viscosity rj. A two-zone mean-field hydrody- 
namics description where the local viscosity is fixed was 
proposed by Zwanzig |l5| ] and in a different context by 
Goodstein to explain ion mobility in superfluid helium 
]l6[ . This picture was elaborated for glasses by Hodg- 
don and Stillinger Much of their calculation can 

now be taken over, but with an additional final averag- 
ing being required. If the viscosity has one value r\i in 
the inner zone and a different value r/ in the outer zone, 
for a incompressible fluid with low Reynolds number, the 
fluid velocity Vi is determined by the linear Navier-Stokes 
equation, 

dp 
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with two values of the viscosity as shown in Fig.(l). p is 
the pressure, and the equation of continuity, 



V-v = 0, 



(6) 



must be satisfied. For a probe particle with radius R 
centered in such a spherical domain with size L, solving 
Eq.(5) and (6) with the proper boundary conditions at 
r = R and r — L, Hodgdon and Stillinger as well as 
Goodstein find, 
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where / is the translational drag force on the probe. This 
is obtained by integrating the pressure over the surface of 
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the sphere, u is the fluid velocity far from the fluctuating 
domain, and C5 is a constant determined by the ratios 
of inner /outer viscosity and L/R. Using slip boundary 
conditions, at the surface of the probe molecule, 



C 5 = l[-3 + 3l 5 + C(3 + 2l 5 )]/d s , 
d s = -2 + 2l 5 + C(4 - 31 + I 5 + 

+( 2 {l-l) 3 (l + l)(2 + l + 2l 2 ) : 



where I — L/R and C = Vo/Vi- 

From Eq.(7), we have the translational diffusion coef- 
ficient in a given domain as 
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In heterogeneous supercooled liquid mosaic, the effective 
translational diffusion coefficient for the whole system 
must be averaged over the distribution of local relaxation 
times, 



(10) 



A = / D t P(s c )ds c . 



The deviation from the homogeneous hydrodynamics re- 
sults, defined as the ratio between "real" diffusion coef- 
ficient and the one predicted by the macroscopic Stokes- 
Einstein Relation, Eq.(I), is 
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With Eq.(3)and (4), we may express £ and I in terms of 



to be expected for the approximated two-zone model we 
applied. It could be connected by a very modest change 
(~ 10%) in the local domain size and doubtlessly reflects 
the crudeness of the simple geometry we have chosen. 
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FIG. 2. The deviation from Stokes-Einstein Relation 
highly depends on the size of the probe particle as shown on 
the graph. The radius of the probe particle is expressed in the 
unit of molecular diameters in the liquid. For a probe with a 
size of several molecular distance, the deviation drops rapidly, 
i.e., the Stokes-Einstein Relation is retained. The points are 
from experiments [I] and the solid line is the prediction made 
by random first order theory for glass transition. 

The temperature dependence of the deviation from the 
macroscopic Stokes-Einstein law is shown in Fig. (3). The 
measurement was conducted for terracene in o-terphenyl 
M. Again, there is a very good agreement between the 
theory and experiment over a large temperature range. 
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where AF and L are the typical values for free energy 
barrier and domain size respectively, which can be ob- 
tained by substituting into Eq.(3) and (4). Using a 
Gaussian distribution of s c with the given and 5s c , we 
readily calculate the effective diffusion coefficient for a 
finite size probe immersed in the heterogeneous mosaic. 

The experimental data are available only for a very 
limited number of materials due to the difficulty of such 
measurements near T g . (At T g the typical relaxation time 
is 1,000-10,000 seconds.) Using the results from our pre- 
vious papers for the distribution of relaxation times, the 
size dependence of the deviation factor for a probe mov- 
ing in o-terphenyl, the most-studied glass-forming ma- 
terial, at T g can be calculated. The result, which does 
not contain any adjustable parameters, is in excellent 
agreement for small probes as shown in Fig. (2). The de- 
viation from the experimental data for large probes is 
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FIG. 3. The temperature dependence of the deviation from 
Stokes-Einstein Relation for terracene diffusing in supercooled 
o-terphenylis shown. The theoretical results (solid line) agrees 
excellently with experimental data (points) [1]. 

The deviation factors depend on the fragility, D, of 
the liquids if they are all measured at the glass tran- 
sition temperature T g . The degree of the heterogene- 
ity, as measured by <5s c /s^, decreases as the liquid gets 
stronger [Q. For very strong glass- forming liquids (like 
SiC>2 with D = 150), one would expect the macroscopic 
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Stokes-Einstein Relation to hold well. The theory pre- 
dicts large deviations only for very fragile liquids (with 
D approximately smaller than 20) as shown in Fig. (4). 
The plot is for self-diffusion (or probe with similar size 
with the liquid molecule). Unfortunately presently, most 
of the experiments have been reported only for a very 
narrow range of substances with D ~ 10. We present the 
fragility dependence thus as a prediction to encourage 
investigations over a wider range of liquids. 
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FIG. 4. The fragility of the supercooled liquids has impor- 
tant effect on the observed Stokes-Einstein Relation deviation. 
For fragile liquids like o-terphenyl, such deviation is as large 
as 3 orders of magnitude, while in strong liquids like Si02, it 
should be only minor. The plot is for self-diffusion at glass 
transition temperature T g . 



proximation we use will hold well. However, for probes 
with very large size comparable to the typical domain 
size, such an approximation may pose a problem as the 
diffusion of rubrene in o-terphenyl indicates. When the 
probe particle occupies multiple domains, an additional 
average over these domains should be taken to get more 
accurate results. Nevertheless, the model developed in 
this paper should serve as a good starting point even for 
larger probe diffusion. 

In summary, the mosaic picture of the supercooled liq- 
uids based on the random first order transition theory 
explains the unusual decoupling of translational diffusion 
from the macroscopic hydrodynamics of glass-forming 
liquids. The picture has as an essential element the 
new length scale of dynamically correlated regions. This 
length is determined to a good approximation by a micro- 
scopic theory of the glassy state itself. The mosaic pic- 
ture also has interesting consequences for the frequency 
dependence and nonlinearity of the response of a translat- 
ing probe. Similarly the mesoscale correlation will affect 
rotational hydrodynamics and relative motion of probes 
in the liquid as is important for diffusion controlled reac- 
tions. We leave these problems for the future. 
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hydrodynamic view of molecular motion in the modern 
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There are a couple of important points to be made clear 
before we conclude our remarks. One is that we use the 
"average" or the measured viscosity as the outer viscos- 
ity rj in the calculation. A more thorough theory should 
treat simultaneously the fluctuations inside and outside 
the domain. Nevertheless the approximation to sepa- 
rate these two terms turns out to be a sound one. The 
physics underlying the simplification is the asymmetry in 
the contribution to the diffusion from the "fast" regions 
and "slow" regions. From the structure of Eq.(10), it is 
clear fast regions with relatively smaller local relaxation 
times (or viscosities) have much larger contributions to 
the averaged effective diffusion coefficient. A little alge- 
bra shows the largest contribution is from regions with 
s c ~ 0.67s^. For these regions the local relaxation time 
is more than 5 orders of magnitude faster than the mea- 
sured relaxation time (about 1,000-10,000 seconds). For 
practical purposes, the outer part of the system acts ho- 
mogeneously for these fast regions. 

The other comment concerns the treatment of smaller 
domains. As we have shown, at T gi the typical domain 
size is about 7 times of the molecular distance. From 
Eq. (4) , we know fast regions are also larger compared to 
the average. For the same reason we discussed above, 
these fast (or larger) regions have a much greater con- 
tribution to the calculated diffusion. The two-zone ap- 



[1] M. T. Cicerone and M. D. Ediger, J. Chem. Phys. 104, 
7210 (1996). 

[2] T. R. Kirkpatrick and P. G. Wolynes, Phys. Rev. A 35, 
3072 (1987). 

[3] T. R. Kirkpatrick, D. Thirumalai, and P. G. Wolynes, 

Phys. Rev. A 40, 1045 (1989). 
[4] X. Xia and P. G. Wolynes, Proc. Natl. Acad. Sci. U. S. 

A. 97, 2990 (2000). 
[5] E. Rossler, Phys. Rev. Lett. 65, 1595 (1990). 
[6] J. Jonas, D. Hasha, and S. G. Huang, J. Phys. C 84, 109 

(1980). 

[7] R. Robinson and R. Stokes, Electrolyte Solutions (But- 

terworths, London, 1965), p. 12. 
[8] F. Fujara, B. Geil, H. Sillescu, and G. Fleischer, Z. Phys. 

B 88, 195 (1992). 
[9] M. Cicerone and M. Ediger, J. Phys. C 97, 10489 (1993). 
[10] M. Cicerone, F. Blackburn, and M. Ediger, J. Chem. 

Phys. 102, 471 (1995). 
[11] J. A. Hodgdon and F. H. Stillinger, Phys. Rev. E 48, 207 
(1993). 

[12] F. H. Stillinger and J. A. Hodgdon, Phys. Rev. E 50, 
2064 (1994). 

[13] M. T. Cicerone, P. A. Wagner, and M. D. Ediger, J. Phys. 
Chem. B 101, 8727 (1997). 



4 



[14] X. Xia and P. G. Wolynes, Cond. Matt. Preprint 0008432 
(2000). 

[15] R. Zwangzig, Chem. Phys. Lett. 164, 639 (1989). 
[16] D. Goodstein, Phys. Rev. B 15, 5210 (1977). 



5 



